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Abstract 

Hierarchical Equal Area isoLatitude Pixelization (HEALPix) system is defined as a layer 
which is used in modern astronomy for using pixelised data, which are required by developed 
detectors, for a fast true and proper scientific comment of calculation algoritms on spherical 
maps. HEALPix system is developed by experiments which are generally used in astronomy 
and called as Cosmic Microwave Background (CMB). HEALPix system reminds fragment of 
a sphere into 12 equal equilateral diamonds. 

In this study aims to make a litetature study for mapping of the world by mathematical 
equations of HEALPix celling/tesellation system in cartographical map making method and 
calculating and examining deformation values by Gauss Fundamental Equations. So, although 
HEALPix system is known as equal area until now, it is shown that this celling system do not 
provide equal area cartographically. 
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1. INTRODUCTION 

Map projection, which is classically defined as portraying of world on a plane, is effected by 
technological developments. Considering latest technological developments, we can offer 
definition of map projection as "3d and time varying space, it is the transforming techniques 
of space's, earth's and all or part of other objects photos in a described coordinate system by 
keeping positions proper and specific features on a plane by using mathematical relation." 

Considering today's technological developments, cartographical projections are developing 
on multi surfaces (polihedron) and Myriahedral projections which is used by Wijk (2008) for 
the first time. Within this scope, HEALPix which is developed for fast, true and scientific 
commentary on global sky mapping and uses unreal calgoritm is proposed for cartographic 
studies. 

In this study Gorski et al. (1998b) and (2005) and HEALPix which is defined as usage of 
multipying the globe to 12 squares, will show that Calabretta's design (2004) and Calabretta 
& Roukema (2007), which is said to area preserving and based on mathematical cell 
pixelisation, is not area preserving. 

2. HEALPix 

HEALPix is a software which has unreal algoritms and visual presentations, which is required 
via developed detectors in modern astronomy and used for scientific commentary of Cosmic 
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Microwave Background (CMB) data on global maps fastly and truely. HEALPix can be 
defined as a sky map format and visualisation analysis mass between information 
technologies and space units and has a wide range of use. 

The best results are get when it provides below features for global sky mapping: 

a- Hierarchical database usage 
b- Equal Area 

c- Isolatitude feature for harmonical assessments. 

A single HEALPix provides above features and its naming is based on these three 
specifications. 

Mathematical requirements which are needed by CMB are the basics of HEALPix (Gorski et 
al., 2005). Mathematical unit sphere [0, 1]*[0, 1] on in the basic of HEALPix is based upon 
mapping of each pixel in K*H. Here, K shows the pixel number between north and south 
poles, H shows pixels aroung the pole. As stated by Gorski et al., (1998b), O'Mullane et al., 
(2000) and Gorski et al., (2005) pixel borders are no geodesic and as pixel corners are not big 
circle arcs, pixels are not squares. 

Considering sky map data standard, Gorski et al., (2005) has chosen H = 4, K = 3 as standard 
in his study in 2005. Total pixel number in standard application is 12. In other words, square 
is divided into 12 equal pixels. Standard choice; 

a. Having less than 4 pixels in poles for avaiding narrow angles. 

b. Fast harmonical applications in equatorial areas. 

c. Minimising pixel shape change in equatorial areas to provide equal area feature. 

HEALPix is understood as 12 devided pixels to give desired resolution and pixels 0-3 north, 
4-7 equatorial and 8-11 south pixel numbered. Standard HEALPix application of globe to 
plane is seen on Figure 1. Calabretta (2004) and Calabretta & Roukema (2007) created this 
projection's cartographical application of such a study. 




0.0 0.5 1.0 1.5 2.0 

x /n 

Figure 1. Showing Globe HEALPix on Plane (Gorski et al., 2005) 
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2.1. Cartographical Application 

Cartographical application of HEALPix, is characteristically defined as below by Calabretta 
(2004) and Calabretta & Roukema (2007). 

a. It is equal area; area on the globe is transferred to projection plane wtihout changing. 

b. Parallels are horizontal plane lines as 3y/d<J> = 

c. Parallels are always divided in the same way. 

d. Pixels are in the area preserving equilateral quadrangles and considerd as ±180° in 
latitude-longitude network. For latitude longitute (<J>, X) and for the cartesian 
coordinates on projection plane (x, y) are used. 

In creating projection equalities, first Roukema & Lew (2004), Calabretta (2004) and 
Calabretta & Roukema (2007) stated the total areas 1/6 north and 1/6 as south pole, and the 
left 4/6 as equatoral area. This situation is HEALPix in Figure 2, and (K, H) = (3, 4) standard 
application. <J>* is transition latitude between equatorial and polar regions. Polar and 
equatorial areas are restricted according to <J>* transition latitude. 



Figure 2. HEALPix (K, H) = (3, 4) applicaiton's standards illustration (Roukema & Lew, 
2004). 

Considering above stated definitions, HEALPix is the combination of Collignon projection's 
(1865) in pole areas, and Equal Area Cylindirc Projection of Lambert (1772) in equatoral 
parts (Collignon, 1865; Tissot, 1881; Lutque & Matarazzo, 2004; Calabretta & Roukema, 
2007). HEALPix projection is not seen as documentation in cartography documents and not 
positioned in a web search. Especially, it does not take place in Snyder's (1993) revision of 
cartography history. MicroCAM web site of Illionis State University, commission 
cartography projections of international cartography union presents 320 projections in a 
catalogue (Anderson, 2003; Calabretta, 2004; Calabretta & Roukema, 2007). None of these 
projections have any resemblance with HEALPix projection. Equations which are created 
from these projections are used and Jacobian determinant and Gauss Fundamental Equations 
are used. 
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Figure3. EdouardCollignonprojections (Calabretta&Roukema, 2007). 
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Figure 4. Lambert's Equal Area Cylindiric Projection (Yerci, 1997) 



Gorski et al. (2005), has chosen K = 3, H — 4 application as standard application of numeric 
sky mapping in 2005. Calabretta (2004) and Calebretta & Roukema (2007) use K = 3 
application in their studies. 



2.2. Transition Latitude 4>* Calculation 



On sphere, steradian angle value which sees should be 2%/3 which is 1/6 of the 4ti 
steradian angle as seen in Figure 2. According to this information calculation is required; 

f 271 d0 C! 2 cos &d& = 2n/3 

^ 27rsin[0]^ 2 = 2?r/3 
=> 2n(l - sintfO = 27r/3 

^> sin0 t = 2 / 3 (1) 
=> sin" 1 2/3 -41.8° (2) 
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2.3. Coordinate Calculations 

HEALPix structure use octahedron as base (O'Mullane, 2000). Calabretta (2004) and 
Calebretta & Roukema (2007) configurated reference point as local coordinate systems origin 
and; 

(PoMhealpix = (0- 0) (3) 

In process, pixel's North, Equatoral and South parts should be checked and reletad parts 
formulas should be used. 

2.3.1. Equatorial Region 

In equatorial region \<P\ < <£* and Lambert Equal Area Cylindirc Projection is used. 
Calabretta (2004) and Calebretta & Roukema (2007) gives relations to be used in equatorial 
region as below. 

x = — sinfp (4) 

y = A (5) 

2.3.2. Polar Regions 

In pole regions (> <J>*)and Collignon projection is used. Calabretta (2004) and Calabretta & 
Roukema (2007) gives relations to be used in pole region as below. 

* = (« 
y = X c + {X-X c ]o (7) 
Here, 

a = Vff(l-|sin<P|) (8) 

a is the ratio of the distance of the pole from the paralel of O to that of the pole from the 
paralel transition latitude O*. 

^ 2 U^1S00)H l-o, + \_ 1BQ0 (9) 
c H \ L 360° 2 J / v ' 

^ _ rl ... K uneven number if <P > 
^0... otherwise 

Numerical results are get from (9) numbered relation. These formulas in South hemisphere 
are applied in equator by mirror symetry. 
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3. APPLICATION 



We will check HEALPix's equal area feature by Gauss Fundamental Equation. At the same 
time we will map HEALPix standard application which is produced by us in VoteMap 
software, we will study on differences between spatial differences of pixels. 

3.1. Control of Equal Area Feature by Gauss Fundamental Equation 

The earth is either considered ellipsoid or sphere, it is impossible to transfer details on it to a 
map plane without deterioration and this deterioration is called as deformation. Deformation 
issues are mentioned by French Tissot (1881) for the first time. 

• h: coefficient of deformation in the direction of meridian cycle (North-South direction) 



• k: coefficient of deformation in the direction of parallel (/c = r{) 

shows and there is two perpendicular direction which length deformations reach extreme 
values (Ucar et al., 2004). 

• a: Maximum length deformation coefficient (a = r max ) 

• b: Minimum length deformation coefficient {b = r min ) 

The directions in which length deformations are maximum and minimum are called as "Main 
Deformation Directions". On original surface, parameter curves, intersect each other 
perpendicularly. In case they intersect perpendicularly on the surface of projection, main 
deformation directions and parameter curve deformations are equal to each other. So; 



Unit disc on the sphere transforms into ellipse shape while being transferred to projection 
plane. So transforming of a specific area to projection plane without being changed is called 
as equal area feature. Equal area condition which takes place in Gauss Fundamental 
Equatitons in Maling (1992) is used as "First Fundamental Quantities"; 




(ID 




(12) 




(13) 



In order to have equal area, it should be; 



hxk = axb = 1 



(14) 



3.1.1 . Equatorial Region 



If (4) and (5) numbered equations are put in (12) and (13) numbered formulas; 




(15) 
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fx = (16) 

£=° < 17 > 

|=1 (18) 
If values required in (12) numbered equation are put in (15) and (17) equations; 

(Ktc \ Ktc 

— COS0) -» h = — cos0 (19) 
2H ) 2H v ' 

is required. If values which are required from (13) relation are put in (16) and (18) relations; 
G = 1 -» k = 1/cos (20) 
If values which are required from (14) relation are put in (19) and (20) relations; 
h x k = 1 

=> — COS0 X = 1 

2/f cos <J> 

^> — = 1 

^H = ^ (21) 

H value is required from (21) numbered relation which provides equal area feature. When K = 
3 is chosen in HEALPix standard application, 

H = y = 4,71238898 (22) 

is found. When H = 4,71238898 and H = 4 is chosen, calculated coordinates are shown in 
Table 1. When H — 4,71238898 is chosen,it is seen in Table 2 that there is no area 
deformation. However, in application it is processed as H = 4 and naturally, area deformation 
occurs. 

3.1.2Polar Regions 

If values which are required from (12) and (13) relation are put in (6) and (7) relations; 

n n 

— = -x , H x cos0 =- x K x cos0 = — cos0 (23) 

d& 2 Vtf(l-sin4>) 2 a 2Ha v ' 

fx = (24) 

5y = _« X (±M XC0S(P (25 ) 
dt> 2 a v ' 
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If (A — l c ) = is accepted to fulfill condition of = 0, (25) numbered relation would be 
as below; 

5y = _^ x O z A £ ) xcos0 = () 

30 2 a ' 

< 27 > 

If values which are required from (23) and (26) relation are put in (12) relation; 

(v — \ 2 — 

- x x cos <Z> J -h>/i = -^-x-^-xcos<Z> (28) 

is required. If values which are required from (24) and (27) relation are put in (13) relation; 
G = o 2 ->k = — (29) 

COS 

When reqiured values from (26) and (27) relations are put in (12) numbered equation by 
condition of area preserving; 

h x k = 1 

Kn cos * cr _ . 
— ^ X — 1 

2Ha cos* 

2H 

=> H = y (30) 

H value which provides equal area is required from (30) numbered relation. When K = 3 is 
chosen in HEALPix standard application; 

H = y = 4,71238898 (31) 

When H = 4,71238898 ve H = 4 is chosen, calculated coordinates are shown in Table 1. 
When H = 4,71238898 is chosen, it is shown in Table 2 that there is no area deformation. 
However, in the application it is applied as H = 4 and naturally deformation occurs. 
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Table 1. Cartesian coordinates (x, y) for geographical coordinates for HEALPix standard 
application (K,H) = (3,4) and (K,H) = (3,4.71238898) 



Geographical 
Coordinates 



<p = 50° 
A= 15° 



H =4 K = 3 




x = 
y = 0,261799387 
x = -0,304914004 
= -0,261799387 
x = 0,304914004 
y = -0,261799387 



x = 0,912809511 
y = 0,346740286 



H = 4,71238898 K = 3 





x = 
y = 0,261799387 
= -0,258819045 
= -0,261799387 



x = 0,258819045 
y = -0,261799387 




x = 0,774816778 
y = 0,250434102 



<P 
A 



-50° 
15° 



x = -0,912809511 
y = 0,346740286 



x = -0,774816778 
y = 0,250434102 



q> = -50° 
A= -15° 



x = -0,912809511 
y = -0,346740286 



x = 1,570796327 



x = -0,774816778 
y = -0,250434102 



3.2 VoteMap 



VoteMap programme is created by C++ and not only for one projection but for all projections 
which can be expressed by mathematical equatitons. In programme, data about the chosen 
projection, and related mathematical equations, parameters, statics and drawing is done by 
standard values which are valid for all projections. 



When VoteMap programme is used in HEALPix (K,H) = (3,4), cartographic ally view of 
Figure 5 is required. By the help of VoteMap programme, as in Calabretta (2004) and 
Calabretta & Roukema (2007) 12 pixels studies, Figure 6 view is obtained. It is seen that in 
Figure 7 as pixels are not exactly equilateral quadrangles, transition latitude, creates a 
trapezoid with 45 degree latitude and this value is called as "Overlap Angle" by us. Using 
Figure 6 and Figure 7, we can see the deformation differences in regards of km 2 in Table 3. 
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Table 2.Cartesian coordinates (x, y) for geographical coordinates for HEALPix 
application^, H) = (3, 4.71238898) and Area Calculation on Sphere 



Geography 
Coordinates 


Cartesian 
Coordinates 
(x,y) 


Area Calculation 
X - " 

ZA — 7 l A n r n+l ~ r n A n+l) 


Area Calculation on Sphere 

, A (Pi - 02 02 + 0\ 

F = 2/? -sin cos 

p 2 2 


1(30°, 15°) 


Xi = 0,5000000 

y-t — \JjLi\Jl-l JjDOl 






2(30°, 30°) 


x 2 = 0,500000 


063141026 


063141026 


3(15°, 30°) 


x 3 = 0,258819045 






4(15°, 15°) 

' 


x 4 = 0,258819045 

y 4 — U,i01/ 77JO/ 




1 ' 










1(40°, 15°) 


x 1 = 0,642787609 

_ n T/r-i 700007 






2(40°, 30°) 


x 2 = 0,642787609 


w, yj ■_> 1 j j. / 1/0 


0373R1 708 


3(30°, 30°) 


x 3 = 0,5000000 

— U,JLJJ70/ /D 






4(30°, 15°) 


x 4 = 0,5000000 
v — 761 799387 














1(75°, 15°) 


Xi = 1,120184954 
y 1 = 0,214140029 






2(75°, 30°) 


x 2 = 1,120184954 
y 2 = 0,297843202 


0,026153869 


0,026153869 


3(60°, 30°) 


x 3 = 0,910683602 
y 3 = 0,4022130281 






4(60°, 15°) 


x 4 = 0,910683602 
y 4 = 0,23615612 








10 




Figure 6. 12 pixels view of HEALPix and its required view from VoteMap 
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Figure 7. Showing of Overlap Angle in HEALPix (K, H) = (3, 4) application 

Table 3. HEALPix (K, H) = (3, 4) Deformation, Polar and Equatorial Pixel Difference 
Values. 



H = 4, K = 3 


a 


16,054619 


b 


1,185253 


a*b (Area Deformation) 


18,849556 


Polar Area 


68445,5 


Equatorial Area 


68040 


Areal Difference 


405,5 


£ Polar Area 


547564 


E Equatorial Area 


272160 


E Area 


819724 


E Areal Difference 


1622 


Ratio 


0,001979 


(E Areal Difference / E Area) 




Overlap Angle 


2,378049 


E Areal Difference (km 2 ) 
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4. RESULTS 



As a result of this study about HEALPix (K, H) = (3, 4) which is mostly used variation in 
astronomical cartography studies and created by mathematical base of software is studied. As 
seen in studies, HEALPix does not provide equal area cartographically. At the same time, it 
does not consist of same shaped and equilateral quadrangle pixels. 

As a result of our studies on HEALPix, we think that HEALPix {K, H) — (3, 4) application 
does not provide desired features and should not be used. As we think that all studies has a 
contribution to scientific development, we believe that it will contribute to cartographical 
applications in future. As a result of our studies we believe that Hierarchical Equal Area 
IsoLatitudes Pixelisation (HEALPix) should be renamed as Hierarchical isoLatitude 
Pixelisation (HiLPix). 

5. SOURCES 

Anderson, P., B., 2003. MicroCAM (Computer Aided Mapping) web site (last updated 2003 
Aug), http://www.ilstu.edu/-misrocam/map projections/ 

Calabretta, Mark. R., 2004. Mapping on the HEALPix grid, Astronomy & Astrophysics. 

Calabretta, Mark. R., Roukema, Boudewijn, 2007. Mapping on the HEALPix grid, Mon. 
Not. R. Astron. Soc. 381, 865-872. 

Collignon, E., 1865. Journal del'Ecole Polytechnique, 24, 125. 

Gorski, K., M., Hivon, E., Banday, A., J., Wandelt, B., D., Hansen, F., K., Reinecke, M., 
Bartelmann, M., 2005. HEALPix: A Framework for High-Resolution Discretization and Fast 
Analysis of Data Distributed on the Sphere, The Astrophysical Journal, 622:759-771, April 1. 

Luque, M., Matarazzo, G., 2004. Projection Collignon a meridiens et paralleles 
restilignes, http://melusine.eu.org/svracu-se/mluque/mappemonde/doc- 
collignon/collignon.html . 

O'Mullane, W., Banday, A., J., Gorski, K., Kunszt, P., Szalay, A., 2000. Splitting the sky- 
HTM and HEALPix, MPA/ESO/MPE Joint Astronomy Conference Mining The Sky. 

Roukema, Boudewijn, F., Lew, Bartosz, 2004. A Solution to the Isolatitude, Equi-area, 
Hierarchical Pixel-Coordinate system, The Astrophysical Journal. 

Snyder, J., P., 1993. Flattening the Earth, Univ. Chicago Press, Chicago. 

Tissot, A., 1881. Memoire sur la representation des surfaces et les projections des cartes 
geographiques, Gauthier-Villars, Paris. 

Ucar, Dogan, Ipbuker Cengizhan, Bildirici, I., Oztug, 2004. Matematiksel Kartografya 
Harita Projeksiyonlari Teorisi ve Uygulamalari, Atlas Yayin Dagitim, istanbul. 



12 



Wijk, Jarke, J. van, 2008. Unfolding the Earth: Myriahedral Projections, The Cartographic 
Journal, Vol. 45, No. 1, pp. 32-42. 

Yerci, Mehmet, 1997. Harita Projeksiyonlari Ders Notlari, Selcuk Universitesi Muhendislik- 
Mimarhk Fakultesi Yayinlari, Konya. 



13 



